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Abstract 

Rapidly rotating spherical kinematic dynamos are computed using the combination 
of a quasi geostrophic (QG) model for the velocity field and a classical spectral 3D 
code for the magnetic field. On one hand, the QG flow is computed in the equatorial 
plane of a sphere and corresponds to Rossby wave instabilities of a geostrophic inter- 
nal shear layer produced by differential rotation. On the other hand, the induction 
equation is computed in the full sphere after a continuation of the QG flow along 
the rotation axis. Differential rotation and Rossby-wave propagation are the key in- 
gredients of the dynamo process which can be interpreted in terms of aft dynamo. 
Taking into account the quasi geostrophy of the velocity field to increase its time 
and space resolution enables us to exhibit numerical dynamos with very low Ekman 
(rapidly rotating) and Prandtl numbers (liquid metals) which are asymptotically 
relevant to model planetary core dynamos. 
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1 Introduction 



The magnetic field of the Earth is produced by a dynamo effect in the metallic 
liquid core of our rotating planet. Many efforts have been made successfully in 
the last decade to describe the mechanism of a self induced magnetic field ei- 
ther with experimental models [1,2,3] or numerical simulations [4,5,6,7]. Both 
approaches have limitations. No experiment has been done in rotation while 
rotation is seen as a key ingredient by geophysicists to explain the geometry 
and the amplitude of the geomagnetic field [8]. All numerical models [9,10] 
have introduced the Coriolis force in solving the Navier-Stokes equation and 
the quasi geostrophy (two dimensionality imposed by the Taylor Proudman 
theorem [11]) of the flow participates in the generation of the magnetic field. 
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Thermal convective vortices aligned with the rotation axis are associated to 
surface patches of magnetic field [7] and spatio-temporal behaviors of mag- 
netic and vorticity field are similar. This effect is a direct consequence of the 
prescribed magnetic Prandtl number (P m = v/rj, where v is the kinematic 
viscosity and r\ the magnetic diffusivity) in the simulations. The current com- 
puter capacities limit the computation to magnetic Prandtl number of the 
order of unity [9] while liquid metals exhibit magnetic Prandtl number lower 
than 10~ 5 , even in the planetary core conditions [12]. 

In this paper, we propose an approach that aims at computing very low mag- 
netic Prandtl number dynamos taking advantage of the quasi-geostrophic be- 
havior of the velocity field. For very low Ekman number (E = u/QR 2 , where 
Q is the rotation rate of the spherical container, and R its radius), a quasi- 
geostrophic (QG) approach models correctly the flow in a rapidly rotating 
sphere [13,14]. It consists of the integration of the flow equations along the ro- 
tation axis. Even if the numerical resolution is done with a stream function in 
the equatorial plane (2D), the top and bottom boundary conditions are taken 
into account through slope and Ekman pumping effects. In the context of 
the study of thermal convection in rapidly rotating spherical shells, Aubert et 
al. [15] have compared successfully their QG results with 3D calculations [16] 
and experimental measurement [17]. Low value of P m may imply a separation 
in term of scales and frequencies, between the velocity and magnetic fields in 
a metallic dynamo. This idea has already been applied to kinematic dynamo 
computations at low P m [18]. 

In this work, we compute the QG flow in the equatorial plane with a fine 
spatio-temporal resolution and the velocity is extrapolated to a coarse 3D 
spherical grid where the induction equation is solved. 

In order to demonstrate the validity of this approach, we have decided to apply 
it to a simple case. Instead of a thermal convective flow for which heat trans- 
port has to be modeled, we model the instabilities of an internal geostrophic 
shear layer. This layer, known as the Stewartson layer, is produced by a dif- 
ferentially rotating inner core in a rotating sphere and consists of two nested 
viscous shear layers [19,20]. For a critical Ro number (Ro = AQ/Q, where 
AQ is the differential rate of rotation of the inner core), the Stewartson layer 
becomes unstable [21] and generates Rossby waves [22]. 

As we will show in this paper, such kind of flow can generate and sustain 
a magnetic field. The QG-model allows us to compute dynamos at very low 
Ekman (down to 1CT 8 ) and Prandtl numbers (as low as 3 1CT 3 ). 
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Fig. 1. Sketch of the split sphere geometry. The differential rotation produces an 
axisymmetric Stewartson I? 1 / 4 shear layer which is cylindrical and aligned with the 
rotation axis z. 

2 The equations 

2. 1 Hydrodynamics 

Let us consider a sphere of radius R filled with an incompressible liquid metal 
of viscosity v and magnetic diffusivity rj. The sphere is rotating at Q along 
the z-axis of a cylindrical reference frame (e s ,e^,e z ). The sphere is split at 
the colatitude ± sin -1 (R s /R) (R s /R is set to 0.35). The two polar caps are 
differentially rotating at AQ as shown in figure 1. is chosen to scale the 
time, R the length, (fx pY^ 2 Rn the magnetic field. For low Ekman and Rossby 
numbers, the flow is quasi geostrophic [11]. Taking the curl of the Navier- 
Stokes equation and averaging along the rotation axis z (noted by an overbar), 
we get the QG equation for the z-component of the vorticity u = e z • V x u, 
provided that u s and are independent of z [22]. 

du dui usduj . ,du z — r . — ,_, . 

a + "-a7 + f^- (2 + ^ = Vx(jxB »' e * + £A ^ (1) 

The Coriolis term needs the evaluation of We deduce that u z is a linear 
function of z from the averaged mass conservation equation. Consequently, 
its vertical derivative may be deduced from the non penetration boundary 
condition (J3 effect) and the viscous boundary condition (the Ekman pumping 
effect) [22]. It gives: 



-^ = E 1 / 2 P(u s ,u^s)+P(s)u s (2) 
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where (3(s) = j ^ _ l and L(s) = y/l — s 2 is the half height of a column of 
fluid and 

P(W S ,W0, S) = — ~~TTTT77 + 



2(1 - s 2 ) 3 / 4 



90 



5s 



2(1 - s 2 ) 3 / 2 



(3) 



is the pumping boundary condition in a rigid sphere deduced from Greenspan's 
formula [22]. 

The axisymmetric flow is computed directly from the velocity equation. 

dJ sr + ("-i?) + ^ + 2 <«■> = <CT^> + E ( A w - ¥) (4 > 



where ( ) stands for the 0-average operator. Rigid boundary conditions are as- 
sumed for the velocity at s = 1. For s < R s /R, the top and bottom azimuthal 
velocity are imposed as = sRo. The velocity field is computed using a gen- 
eralised stream function in the equatorial plane as in [22] which guarantees 
the 3D mass conservation. The stream function is expanded in Fourier compo- 
nents along the <fi component. It may be interesting to introduce the Reynolds 
number Re = RoE" 1 directly related to the two controlling dimensionless 
numbers E and Ro. 

In this paper, as a first step, we will only consider kinematic dynamos and the 
magnetic terms in (1, 4) will be neglected. 



2.2 Induction equation 



The velocity field computed with equations (1, 4) in the equatorial plane is 
extrapolated to a spherical grid (on Gauss collocation points) in the physical 
space. This is a straightforward process because u s and are independent 
of z and u z is a linear function of z. Then, the velocity field is changed into 
spherical coordinates (e r , e#, e^) to compute the non linear induction term. 
The dimensionless equation of evolution of the magnetic field is : 

flD 

— = V x (u x B) + P-'EAB (5) 



Changes of magnetic Prandtl number P^ 1 would change directly the magnetic 
Reynolds number R m = ReP m = RoE~ 1 P m which is more commonly used in 
dynamo modeling. This equation is solved using spherical harmonics where 
the magnetic boundary conditions are easy to write [8]. The induction part 
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Fig. 2. z-vorticity maps in the equatorial plane, (a) and (b) : E = 10~ 6 , Ro = 0.0096 
and Ro = —0.0111 respectively. It shows the flow at the onset of hydrodynamic 
instabilities for both signs of the Rossby number, (c) and (d) : E = 10~ 8 , Ro = 0.02 
and Ro = —0.02 respectively. It shows a typical view of the "turbulent" regime for 
Rossby numbers about 30 times critical. The color bar gives the local vorticity scale 
for (c) and (d) only. 

of the code has been checked using kinematic dynamo results [23] and the 
dynamo benchmark [24]. 

2.3 Numerical implementation 

A finite difference scheme is used on an irregular radial grid (denser in the 
Stewartson and Ekman layers). A semi implicit Crank-Nicholson scheme is 
used for linear terms in time whereas an Adams-Bashforth procedure is imple- 
mented for non linear terms. For low P m , cylindrical and spherical radial grid 
steps may differ by a factor 20. Similarly, time steps for the induction equation 
may be much longer than the velocity time steps (as much as 20 times). For a 
run at E = 10~ 8 , the stream function is computed on a cylindrical mesh made 
of 600 radial points and expanded in Fourier series up to degree m = 170 while 
the magnetic field is expanded in spherical harmonics (L max = 79, M max = 32) 
with an irregular radial grid of 150 points for P m = 10 -2 ' 5 . By increasing M max 
and reducing the time step factor, we checked that such truncatures do not 
influence the onset of dynamo action. 



3 Hydrodynamics 

For low Rossby numbers, the split at the spherical boundary produces an 
internal shear layer in the fluid on a cylinder of radius R s aligned with the 
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Fig. 3. Spectra of the kinetic energy, and both toroidal and poloidal magnetic en- 
ergy for E = 10 -8 , Ro = 0.02 (30 times critical) and P m = 5 10~ 3 (equivalent to 
Re = 210 6 ). 



rotation axis. This geostrophic viscous layer consists of two nested layers of 
different widths as revealed by the asymptotic study of Stewartson [19] and 
illustrated later by a numerical study of Dormy et al. [20] ; a larger one of size 
E 1 ^ which accommodates the jump in azimuthal velocity geostrophically and 
a narrower one of size E 1 ^ 3 , ageostrophic, which corresponds to an axial jet 
insuring the mass conservation. 



In our previous study [22], we presented the QG model, which can reproduce 
only the E lj/4: layer, and we studied the linear perturbations of this geostrophic 
internal viscous layer. It becomes unstable when the Rossby number exceeds 
a critical value Ro c which varies as (3E 1 / 2 [22]. At the onset, the instability is 
a Rossby wave, azimuthal necklace of cyclones and anticyclones of size 

£1/4 

which propagates in the prograde direction as shown in figure 2ab. Super 
rotation (Ro > 0) generates a spiraling flow outside the shear layer while the 
flow is mainly located inside the shear layer for Ro < 0. For supercritical Ro, 
the flow exhibits larger vortices (fig. 2cd) which are time dependent but still 
drifting as Rossby waves. The flow stays mainly concentrated in the shear 
layer. Figure 3 shows the kinetic energy spectra E(k) of this QG turbulent 
flow. It is very steep : E(k) ~ k~ 5 which is the spectrum predicted by Rhines 
[25] for turbulence in presence of Rossby waves. This steep spectrum suggests 
that the small scales of the flow may be neglected in the induction equation. 
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4 Dynamo action 

For a given Ekman number (E = 1CT 6 to 1CT 8 ), we vary the Rossby number 
Ro from critical to a few times critical and we find the critical magnetic 
Prandtl number P m of the onset of dynamo action by trial and error tests. 
As the flow is time dependent, we detect dynamo criticality on long term 
time variations of the magnetic energy. Unlike most of the kinematic dynamo 
models [26], a critical magnetic Prandtl number was found for every set of 
dimensionless numbers (E, Ro) we have computed. In figure 4, we plot the 
calculated critical magnetic Prandtl number as function of the Reynolds 
number Re = Ro/E. As expected, we found that an increase of the forcing 
(Ro) for a given E reduces the critical magnetic Prandtl number. A decrease of 
the critical magnetic Prandtl number is also observed as we lower the Ekman 
number. These two effects may be summarised by the use of the magnetic 
Reynolds number R m . The data points in figure 4 are roughly compatible 
with the line R m = 10 4 . A critical magnetic Reynolds number R c m of 10 4 is 
only indicative because the details of the flow generate large deviations (factor 
3) from this simple law. The minimum critical magnetic Prandtl number of 
0.003 has been found for E = 10~ s and Ro = 0.02. 

The critical magnetic Prandtl number is not independent of the sign of the 
differential rotation (sign of Ro). This is expected because the flow is quite 
different in the two cases as shown in figure 2. A negative differential rotation 
seems to lead to slightly lower dynamo thresholds. 

Antisymmetric axial velocities (u z (z) = —u z (—z)) and symmetric orthoaxial 
velocities (u s ^(z) = u s ^(—z)) generate two independent families of growing 
magnetic field in kinematic dynamos known as the dipole and quadrupole fam- 
ilies [27]. The geometry of the two families are shown in figure 5a and 5b : the 
dipole family is dominated by an axial dipole, whereas the quadrupole family 
exhibit a strong axial quadrupole. Each family has a different critical magnetic 
Prandtl number. As shown in figure 4, we found that the dipole family has al- 
ways a larger critical magnetic Reynolds number than the quadrupole family. 
This result is quite different from the conclusion of the work of Sarson and 
Busse [28]. Using Kumar and Roberts kinematic dynamos, they found that 
prograde spiraling of columns and prograde zonal flows favor dipole magnetic 
fields. 

In both families, the strongest magnetic fields are produced in the Stewart- 
son shear layer deep inside the sphere. The typical spectra given in figure 3 
show that the computed magnetic fields are dominated by both toroidal and 
axisymmetric components. At the surface of the sphere (figure 5c), the radial 
magnetic field is also mostly axisymmetric, and the non-axisymmetric part 
is clearly associated to the geostrophic vortices produced in the Stewartson 
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Fig. 4. Dynamo onset for different parameters : Critical magnetic Prandtl number 
versus the absolute Reynolds number Re = \Ro\E~ 1 . Dipole and quadrupole 
thresholds are respectively denoted by circles and squares while solid and open sym- 
bols represent positive and negative differential rotation. All the points lie around 
the R m = 10 4 line. 

shear layer. 

The geometry of the magnetic field may be understood in term of aVt effects 
[27,8]. A very large toroidal magnetic field compatible with the azimuthal flow 
is converted to a poloidal magnetic field by the columnar flow through an a 
effect. Any non azimuthal component of the magnetic field is transformed into 
an azimuthal component by the strong differential rotation in the Stewartson 
layer by f2 effect. 

If we consider the magnetic Reynolds number associated to the shear flow 
R% = RoP m E~ l and the magnetic Reynolds number based on the vertical 
velocity = u z P m E~ x k~ x where u z is deduced from the calculation, as well 
as the Rhines wave number k [25]. We may compute the dynamo number Dy = 
\J RmRm [27]. Figure 6 shows that Dy stays roughly constant (between 200 and 
300) as the flow becomes more and more vigorous for the quadrupole family. 
The dipole family seems more easy to excite for negative Rossby number. 
As proposed by Robert [27], this feature may indicate that we have aVl' < 
(where £Y is the radial derivative of f2) in the northern hemisphere for Ro < 0. 

Both the geometry and the onset in term of dynamo number Dy indicate 
that we may look the Stewartson QG dynamo as an afl dynamo where the Q 
effect is produced by the azimuthal shear layer and the a effect by the vortex 
necklace. 
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Fig. 5. Growing magnetic field in kinematic dynamos for E = 10 -8 . (a) and (b) 
are meridian cut of the sphere showing the axisymmetric part of the magnetic field. 
The solid lines are the poloidal field lines and the color map represent the azimuthal 
field, (a) shows a quadrupole field obtained at Ro = 0.02 and P m = 0.005. (b) shows 
a dipole field obtained at Ro = —0.02 and P m = 0.003. (c) is a spherical map of 
the radial magnetic field at the surface of the core, corresponding to case (b), the 
dashed line being the rotation axis. The corresponding vorticity fields are given in 
figure 2cd. 
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i . The critical dynamo 



Busse [29] suggested that the Ekman pumping is important for the dynamo 
process : although the /3-effect produces axial velocities, they are out of phase 
with the axial vorticity at the onset of thermal convection in a rapidly rotating 
annulus and cannot contribute to the mean helicity, whereas axial velocities 
due to Ekman pumping are in phase with the axial vorticity. However, the 
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Table 1 

Table of the critical magnetic Prandtl numbers for the different calculations. 



Ekman pumping flow is of order E 1 / 2 , so that the dynamo process proposed by 
Busse becomes very weak when lowering the Ekman number. In addition, when 
artificially removing the Ekman pumping flow in our dynamo simulations, 
we still observe dynamo action with nearly the same threshold. It seems to 
indicate that the /3-effect alone may produce an efficient a effect, without 
requiring an Ekman pumping flow. 

Furthermore, we have not been able to find a critical magnetic Prandtl num- 
ber with a steady flow (either a time averaged flow or a flow with its time 
evolution stopped at a given time). It implies that the time evolution of the 
flow is a key ingredient for dynamo action in these quasi-geostrophic dynamos. 
The propagation of the Rossby waves is required to put in phase the non ax- 
isymmetric magnetic fields and velocities in order to produce a axisymmetric 
poloidal magnetic field. This type of a effect was proposed in the model of 
Braginsky [30]. Currently, many dynamo experiments are designed with the 
help of numerical simulations (kinematic dynamos). Even if the flow is highly 
turbulent (Re > 10 6 ), mean flow approaches are used for simplicity purposes 
to find the dynamo onset [1,31,32,33]. This method would fail in the case of 
Stewartson dynamos. 
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4-1 Oscillating solution 



As in many aQ dynamo [27], we sometimes obtain a time oscillating solution 
for the Stewartson dynamo. Dipole solutions for E = 1CT 6 do exhibit such 
a behavior. The growth rate of one of these dynamos is plotted on figure 7, 
showing three time scales : the smallest one is the time scale of the velocity 
fluctuations. The intermediate time scale is the time needed for the growth 
rate to go from its minimum value to its maximum. The large time scale is 
the period of oscillation, not linked to any time scale of the flow. 

In the context of kinematic dynamos, this behavior correspond to a complex 
eigenvalue of the linear set of equation [27,8,23]. We use the same formalism 
to explain our result. It may be explained by two coupled magnetic modes B\ 
and B 2 . Let assume that the induction equation may be approximated by the 
following system : 



dBi 

~dT 

dB 2 

dt 



■ X 1 B 1 + K 12 B 2 
X 2 B 2 — K 2 i B\ 



(6) 
(7) 



with all real values. For low coupling (Ki 2 K 21 < (Ai — A 2 ) 2 /2) the eigenvalues 
of this system are real, so that the growing solution is a combination of B\ 
and B 2 . This is the case for the quadrupole family at Ro < 0. However, when 
the coupling K 12 K 2 i is sufficiently strong, the eigenvalues are complex and as 
a result the growing magnetic field oscillates periodically between B x and B 2 . 
The intermediate time scale (time for the growth rate to go from its minimum 
to its maximum) is very close to the phase shift between B\ and B 2 , and one 
of the two modes is dominant near the minimum of the growth rate cycle, 
while the other one is dominant near the maximum, with growth rate close to 
Ai and \ 2 . 



The reversal process at work in our simulations is a smooth periodic evolution 
of the magnetic field, but at the surface it appears to be a sudden sign reversal. 
In fact, a reversed poloidal magnetic field is slowly growing inside the Stewart- 
son layer, moving away the initial poloidal magnetic field until it reaches the 
outer boundary. Then, the reversed dipole magnetic field suddenly appears at 
the surface and ultimately the poloidal field reverses at the center. During the 
time oscillation, the axisymmetric toroidal magnetic field patches in the Stew- 
artson layer migrate toward the equator as reversed polarity toroidal fields are 
formed at higher latitudes. This migration could be understood in terms of 
Parker dynamo waves [34,27]. 
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Fig. 7. Magnetic field reversal observed at E = 10 , Ro = —0.08 and P m = 0.1. The 
graph shows the evolution of the growth rate of the magnetic energy as a function 
of time (in magnetic diffusion time units). 
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5 Conclusion 



In summary, we have computed a quasi-geostrophic dynamo based on a Stew- 
artson shear layer flow. The scale separation approach works because the small 
scales of the flow in our rotating sphere are negligible (very steep kinetic en- 
ergy spectrum E(k) ~ A; -5 ). Our preliminary results may be interpreted in 
terms of afl dynamo. The Q effect is done by the shear of the Stewartson layer 
itself whereas the a effect is produced by vortices associated with the Rossby 
waves due to the instability of the shear layer. These understandings are very 
encouraging for our on-going experimental modeling of the geodynamo. As 
described in Car din et al. [35], we are building a spherical Couette experiment 
using liquid sodium which may validate and enlarge our present numerical 
findings. 

For the first time, we have computed a spherical dynamo with a very low 
magnetic Prandtl number (< 1CT 2 ) and a very low Ekman number (1CT 8 ) 
(corresponding to a very high Reynolds number Re > 10 6 ). Even if our di- 
mensionless parameters stay far away from parameters of planetary cores, our 
calculations use dimensionless numbers which are in the correct asymptotic 
regime for the modeling of the geodynamo. The key ingredients of our ap- 
proach is to take into account a specific property of the rotating fluid (QG) 
which allows us to use a 2D model to compute the flow evolution, and the 
separation of scales between the magnetic field and the velocity field, allowing 
us to use a coarse 3D mesh for the magnetic field. 

We also showed that in the case studied in this paper, the mean flow or the 
static flow fails to produce a dynamo while the fully resolved time-dependent 
flow succeeds. Indeed, the time evolution of the flow and the f3 effect are key 
ingredients for dynamo action in our models, while the Ekman pumping can 
be neglected without losing the dynamo effect. 

The next step will be to add the Lorentz force in the QG equation to compute 
saturated dynamos. One of the difficulty is to compute the action of the large 
magnetic field on the small scale motions of the fluid. Preliminary results 
are encouraging and exhibit saturated dynamos very close to the kinematic 
dynamos described here. 

A quasi geostrophic approach could also be used to build thermal convective 
dynamos. A zonal geostrophic flow is produced by the Reynolds stress of the 
thermal columns [17,36] but its amplitude is much lower compared to the 
differential rotation imposed in the Stewartson problem. Would it be enough 
to start an Stewartson dynamo type? for what forcing? Would it work for very 
low Ekman and magnetic Prandtl numbers? 
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